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1. Introduction 



Let us recall that M = 2 supersymmetric (SUSY) extensions of the KdV equation 
have been largely studied in the past [THSj in terms of integrability conditions and 
solutions. Such extensions are given as a one-parameter (a G R) family of Grassmann- 
valued partial differential equations with one dependent variable A(x,t; 61,62) where 
the independent variables are given as a set of even (commuting) space x and time t 
variables and a set of odd (anticommuting) variables 61, 62. The dependent variable A 
is supposed to be an even superfield and satisfies 

A t = {-A xx + (a + 2)AD 1 D 2 A + (a - 1)(D 1 A)(D 2 A) + aA 3 ) xl (1) 

where D\, D2 are the covariant superderivatives given by 

A = 6id x + d 6i , 1 = 1,2, (2) 

and such that D\ = D\ = d x . 

Since the odd variables satisfy 6\ = 6\ = 6162 + #2^1 = 0, the dependent variable A 
admits the following Taylor expansion 

A(x, t; 6 U 6 2 ) = u(x, t) + d^iix, t) + 6 2 &(x, t) + e 1 6 2 v(x 1 t), (3) 

where u and v are commuting (bosonic) complex valued functions and £1 and 6 
are anticommuting (fermionic) complex valued functions. Equation (pD) can thus be 
decomposed set of two bosonic equations 

ut + (u xx - au 3 + (a + 2)uv - (a - l)£i&)x = 0, (4) 
Vt + {Vxx + 3v 2 - 3au 2 v - (a + 2)uu xx - (a - l)u 2 x ) x 

-(366;x + (a + 2)66;,. - 6a«66)x = (5) 

and two fermionic ones 

6;t + (Ci-,xx + 3(w - au 2 )6 - (a + 2)u£ 2i x - (a - l)itx6)x = 0, (6) 
6;t + faycx + 3(v - an 2 )6 + (a + 2)<i ; x + (a - 1)^6)^ = 0. (7) 
We easily see that this last set can be decoupled. Indeed, taking p± = 6 ± ^6; we g e t 

P±-,t + (p± ; xx + 3(u - aw 2 )p ± ± i(a + 2)up±. x ± i(a - l)w x .p ± ) x = 0. (8) 

Solutions of equation (JT|) have been obtained using different methods. Among them, 
an original approach [2] has adapted the classical symmetry reduction method to the 
SUSY context. Indeed in this approach, starting with the SUSY equation ([!]), it has 
been reduced to a set of only one bosonic and one fermionic equations using invariance 
properties associated with invariant superalgebras. 

It looks like a reduction from M = 2 to M = 1 SUSY but it is not quite the same as 
we will see later. It also gives different results than a reduction method which consists 
of setting to zero some of the dependent variables (see [5] for example). 



Let us mention that the bosonization approach [6] has also been used to solve 
such systems. It consists in expanding the bosonic and fermionic fields in n-fermionic 
parameter space. In the following developments, we will consider the simplest case 
(one-fermionic parameter bosonization) which consists of writing the fermionic fields as 
d = C/»(* = 1)2), with fi a bosonic complex function and £ an odd parameter such 
that ( 2 = 0. In the decomposed original system of equations (j3J) and (jSJ), the fermionic 
contributions thus disappear and we get purely bosonic equations: 

u t + {u xx - au 3 + (a + 2)uv) x = 0, (9) 
Vt + {v X x + 3i> 2 — (a + 2)uu xx — (a — l)u x — 3au 2 v) x = 0. (10) 

It is interesting to see that for v = ±iu x , (fTOl is a consequence of (jHj). 

This is not the case for more general bosonization procedure. For example, in the 
2- fermionic parameter space spanned by (j and the bosonic fields u and v takes the 
forms Uq + QiC, 2 ui 2 an d v + CiC2^i2, respectively. We see that the functions u and t> 
satisfy the same set of equations (Q and (fTOl) . On the other hand, the fermionic fields 
£1 and £2 have the new forms dpi + (2P2 and Ci^i + C292, so that the set of fermionic 
equations (jSJ) and ([7j) is doubled. The evolution equations on u± 2 and Vi% involve all the 
bosonic and fermionic dependent variables. 

In the following, we will take a = —2 in equation ([T|) where it is well-known that 
we get solitons solutions as travelling wave solutions prH] and also rational similarity 
solutions as in the classical case 0IE]- Such solutions will be generalized in the SUSY 
case using the symmetry reduction method and also the Hirota formalism. 

Let us write the set of equations we are working with in this case. The equation 
([T]) becomes: 

A t = (-A xx -3(D 1 A)(D 2 A)-2A%, (11) 
and admits the decomposition 

ut + {u xx + 2m 3 + 3^ 2 ) x = 0, (12) 
v t + (v xx + 3v 2 + 6u 2 v + 3u 2 x ) x - (366;, + = 0, (13) 

£i;t + (£i ; xx + 3(v + 2u 2 )^ + 3u x &) x = 0, (14) 
6;t + (6;xx + 3(v + 2m% - 3^6), = 0. (15) 

The equations (JSJ become 

P± ; t + (P±;xx + 3(f + 2n 2 =f iu x )p±) x = 0, (16) 
and thus equations (fT2|) and (fl3l can be written as: 

3 

«t + (u xx + 2u 3 + -ip + p_) x = 0, (17) 

3 

v t + (v xx + 3v 2 + Qu 2 v + 3u 2 x ) x + (-(p_ - p + )(/°-;x - P+;x) - 6iup + p_) x = 0. (18) 



In particular, for one-fermionic parameter bosonization, we get the following system of 
nonlinear partial differential equations : 

u t + {u xx + 2m 3 ),. = 0, (19) 
v t + (v xx + 3v 2 + 6u 2 v + 3u 2 x ) x = 0. (20) 

We see that u satisfies the mKdV equation while v solves a modified KdV equation. 
2. Symmetry reduction method and solutions 

Let us here briefly recall the symmetry reduction method as it has been adapted to 
the SUSY context [2]. The Lie superalgebra of symmetries of equation ([T|) is a (3|2)- 
dimensional superalgebra with three even generators associated with time and space 



translations and dilations and two odd generators. We get explicitly; 

V t = dt, V x = d x ,V = xd x + 3td t + ^(M* + 2 d d2 ) - Ad A , (21) 

Qi = 0id x -d ei , Q 2 = 6 2 d x -d d2 . (22) 
The non zero structure relations are given by 

[V t ,V]=3V t , [P X ,V]=V X , [Qi,V] = ^Qi, z = l,2, (23) 

{QuQi} = {Q 2 ,Q 2 } = -2V X . (24) 



Two distinct sub-superalgebras are considered and used to reduce the equation (QQ) 
to a set of nontrivial ordinary differential equations. The first sub-superalgebra G\ 
is generated by {Vt — cV x , Q± + 1Q2} (c is a real parameter) and corresponds to 
translational invariance that will give rise to the so-called travelling wave solutions. 
The second sub-superalgebra G 2 is generated by {V, Q 1 + iQ 2 } and corresponds to 
dilation invariance that will give rise to the so-called similarity solutions. 

In order to perform the symmetry reduction, we have to find the invariant of 
the action of the corresponding subgroup Q\ or Q 2 on the independent and dependent 
variables and rewrite the equation (CQ) in terms of these invariants. For the subalgebra 
G\ and corresponding subgroup Q\, we get the invariants [2] 

y = x + d + z0 x 2 , 6 = 0i + i9 2 , A = A{y, 6) (25) 

and the reduced equation (after integration with respect to y) is 
A yy - aA 3 + i(a + 2)AA y + cA 

- 2i(a + 2)6AA ye - 2%{a - l)9A e A y + c x + Bk x = 0, (26) 

where C\ and k\ are complex integration constants. It is not a M = 1 SUSY reduction 
since equation (1261) cannot be expressed in terms of the superderivative Dq = 9d y + dg. 



Using the decomposition of the invariant superfield A(y, 9) = u(y) + 9p(y), we get 
the following system of ordinary differential equations: 

u yy — au 3 + i(a + 2)uu y + cu + c\ = 0, (27) 
Pyy ~ "i( a + ^) u Py + ( c — 3cm 2 + i(4 — a)u y )p + k\ — 0. (28) 

Identifying A(y, 0) with t; 6*i, #2) given by equation <^ , for y — x + ct, we get 

u(x,t) = u(y ),v(x,t) = i ^-^- ,Ci(x,t) = -i£ 2 (x,t) = p(y )- (29) 

dyo 

We see that the two fermionic functions reduce to one independent function p since 
p + = and p_ = 2£i = 2p. 

For the subalgebra G 2 and corresponding subgroup £2, we get the invariants 

z = H( x + ie x e 2 ), e = e 1 + id 2 , w = &A(z, e) (30) 

and the reduced equation (after integration with respect to z) is 
W zz - aW 3 - -zW + i(a + 2)WW Z 

- 2%{a + 2)eWW zB - 2%{a - l)6W e W z + c 2 + 9k 2 = 0, (31) 

where c 2 and k 2 are complex integration constants. 

Since the invariant superfield W can be decomposed as W(z, 6) = w(z) + 8(f)(z), we 
get, from ( 13TT) . the following system of ordinary differential equations: 

w zz — aw 3 + i(a + 2)ww z — -zw + c 2 = 0, (32) 

4>zz — i{a> + ^)w(j) z + (— -z — 3aw 2 + i{A — a)w z )(p + k 2 = 0. (33) 

The only difference between the sets of reduced equations ( 1271) and ( 1281) is that the 
constant parameter c has been replaced by —\z in (|32|) and (|3~3"I) . 

Identifying 6 1 ) = t~3VF(2:,#) with v4(x, t; 81, 6 2 ) given by equation ([3]), for 

£ = t 3x, we get 

u(x,t)=t 3 W (zo),v(x,t) = it 3 , - ,£i(x,t) = -i£ 2 {x,t) = t s <f>{z ). (34) 

azo 

As we see both reductions apply to equation (jTJ independently of the parameter a. 

2.1. Soliton travelling wave solution for a = —2 

From f )27j) and ( 128]) . we get the decoupled system of equations as 

Uyy + 2u 3 + cu + C\ = 0, (35) 
Pra + ( c + 6m 2 + 6iu y )p + ki = 0. (36) 



General solutions of this set have been given in [2]. Here we give some particular 
solutions, we will start with to make the connection with the Hirota formalism. 
The mKdV equation ( 1351) admits the particular soliton solution 

u(y) = sech(j/), (37) 

if we take c\ = and c = — 1. The fermionic equation (|36|) thus admits the particular 
solution (ki = and £ is an odd parameter) 

P(y) = (po(y) = C sech 3 (y)(sinhy + -(1 - sinh 2 y)). (38) 

Using fl29|) . the original solution t; 6*i, # 2 ) is thus given by (j3J) with 

u(x,t) = sech(x — £), v(x,t) = —i sinh(x — t) sech 2 (x — t), (39) 

^(x, t) = -i£ 2 (x, t) = C sech 3 (x - t)(sinh(x - t) + ^(1 - sinh 2 (x - t)). (40) 

We see that we have indeed v = iu x . 

Let us observe that using this symmetry reduction, we have in fact solved our 
original problem (fTTj) with the following constraints £i = — 1^2 and v = iu x . This means 
that we have solved the set of decomposed equations: 

u t + {u xx + 2u 3 ) x . = 0, (41) 

ti ; t + (6;xx + 6(iu x + u 2 )^) x = 0. (42) 

Figure 1 shows the behaviour of u as given in f|3T|) and the real and imaginary parts 
of po given in (138]) as a function of y = ( x ~ t). 
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Figure 1. The functions u(yo), Re(po(yo)) and Im(po(yo)). 



Rational similarity solution for the case a = —2 
From (1321 and ( 1331) . we get the following system: 



w zz + 2w 3 - -zw + c 2 = 0, (43) 
o 

+ (-\z + 6w 2 + 6iw z )<p + k 2 = 0. (44) 



The equation (1431) is known as the second Painleve equation and rational solutions for 
the bosonic mKdV equation are known [THTTj. A particular solution is (with c 2 = 

z 3 -6 . d , ,z 3 + 12, 

With such a solution, the fermionic equation (144]) becomes 

36 + z 3 



0, (46) 



when k 2 = 0. Using the following change of variables 

^W=(^(x(z)), xW = 37! zl ' (47) 

such that = 0, we see that ^ satisfies the modified Bessel equation, 

xVxx + X^x " (X 2 + ^ = 0, (48) 
with i/ = ±|. We get, in particular, 

t(z) = CVZl- i (Jj=z*), (49) 

where I^(x) is a modified Bessel function. The original solution A(x, t; 6\, 9 2 ) is thus 
given by equation where 

, . x 3 -6t . . x 6 - 48tx 3 - 72t 2 , . 

U X ' = 2l i 3 I nrt ' U X ' t = 2 2f 3 i 10^2 ' 50 

a;(ar + 12t) ar(ar + iity 

£i(x,t) = -i&(x,t) = H</>(H X ). (51) 

Here again, we see that v = iu x . As before, u and £i solve ( l4Tj) and (1421) respectively. 

Figure 2 shows the behaviour of the imaginary part of u as a function of x and t. 
Such solution of the M = 2 SUSY equation has not been considered before and may 
be compared with the two-soliton one that will be given in the next section using the 
Hirota formalism. 

We refer to Appendix 1 for a discussion of the series of rational similarity solutions 
of equation (1430 constructed using the Yablonskii-Vorob'ev polynomials 0,(11] and the 
corresponding resolution of the fermionic equation 



3. Hirota formalism and solutions for the case a = — 2 

The Hirota formalism is a well known process in the classical and in the M = 1 SUSY 
cases [12HT5]. This formalism has been used, in particular, to obtain iV-soliton solutions. 
We examine the possible generalizations in the M = 2 SUSY case [HE] following the 
idea given in [16]. We write the superfield ([3]) in the following form 



A(x, t- 9i, 6 2 ) = u\x, t- Or) + 9 2 e f (x, t; 9 1 ), 



(52) 



Figure 2. The function Im(u(x,t)) for t = 



-10,0,10. 



where u b , are even and odd superfields, respectively. In fact, comparing with fl3]) , we 
have u b = u + and ^ = £ 2 — Oiv.We thus introduce the following change of variable 
A = d x B, where 

B(x, t; e u 6 2 ) = U\x, t; Ox) + 6 2 rj f (x, t; 9J, (53) 

to equalize the order of the equation with the number of appearence of the x derivative 
in the nonlinear terms. Note that since Df = d x (i = 1,2), Di can be viewed as the 
square root of d x . Equation (11) thus becomes, after integrating once, 

B t = -B xxx - 3(D 1 B X )(D 2 B X ) - 2B 3 X , (54) 

where the constant of integration is set to zero. Inserting the explicit form ( )53|) of B in 
equation ( )5l"l) . we get a set of coupled Af = 1 SUSY equations on U b and rjf: 

U\ + U b xxx - HDiUl + 2(U h x ) 3 = 0, (55) 
4 + VL - HDi4 - m h xx D x U h x + Q^M? = 0. (56) 

Now the strategy is to introduce a change of the dependent variables, U h and rj* , 
in such a way that equations (1551) and (1561) become quadratic. 

Since equation fl55l) is a modified mKdV equation, we use the change of variables: 

Z/ = alog(^), 7^ = /3D 1 log(^) J (57) 

where T\ = t±(x, t; 6i) and r 2 = r 2 (x, t; 6i) are both bosonic superfields and the constants 
a and (3 satisfy a 2 = — 1 and 1 = 1. With no lost of generality, we choose a = — i and 
(3 = 1. 

Equation (IBTI) represents a constraint on the superfield A since we get 

rj f = iD x U h (58) 

and thus 

£f = iD lU b . (59) 
Using fl59l) and comparing ([3]) with fl52l) . we get 

A(x,t;e 1 ,9 2 ) = u(x,t) + (9 1 + i6 2 )£ 1 (x,t) - id 1 6 2 u x (x,t). (60) 



The constraints on the field £2 and v are very similar to the ones obtained in the 
symmetry reduction. We have already noticed that v = ±iu x solves our problem when 
u solves mKdV given in ( 1191) . Here, we see that these constraints are imposed in order 
to produce a bilinear form. 

Inserting equation ( 158]) in equations ( 155]) and ( 156]) . we get 

U b + U b xxx + 2(U b ) 3 = (61) 

and 

D l (U b +U b xxx + 2(U b x f)=0. (62) 

This last equation is trivially satisfied whenever U b is a solution of the mKdV equation 
(16"Tj) . Let us mention that ( 16T|) is still a SUSY equation and if U b takes the decomposition 
Uq(x, t) + 9i(p(x, t), we get the classical mKdV equation for Uq and cp satisfies 

<Pt + fxxx + 6(Uq. x ) 2 (p x = 0, (63) 

for which a particular solution is ip = ( Uq. x where C is an odd constant. Such a result 
is common when we deal with M = 1 SUSY KdV and mKdV equations. 
We thus have a direct bilinearization of equations ( 155]) and (156]) 

(V t + Vl)(r v r 2 ) = 0, 
SV x (r v r 2 ) = 0, 

where 

SV x (t v r 2 ) = (D 6l - D e2 )(d Xl - d X2 ) n n(xi] 6i)r 2 (x 2 ; 6 2 )| a; = :Cl = :C2 ,0 1 =e 1 =e 2 , (66) 

is the super Hirota derivative and D® i = d® i + Qid Xi , i = 1,2. Equations (J&U and (]65j) 
are a natural generalization to the susy case of the classical bilinear form of the mKdV 
equation. 

Let us finally mention that, since we have 

u b = u + 0x6 = d x U b , (67) 

such a bilinearization gives the following set of decomposed equations (together with 
v = -iu x and £ 2 = i£i (p+ = 0)): 

u t + (u xx + 2n 3 ) x . = 0, (68) 

£l;t + (6;xx + 6M 2 ei)x = 0, (69) 

which clearly solve (TTTj) and such that a particular solution of the fermionic equation is 
£1 = C u x- The comparison of this set of equations with the set (HT]) and (H2]) obtained 
in the symmetry reduction, suggest that we will get different solutions for the fermionic 
components £1 and £ 2 of the superfield A. 




3.1. Super soliton solutions 

The Hirota formalism helps us to recover the travelling wave solution or one super soliton 
solution but we also get the N super soliton solutions (HUE]. Indeed, we take 

n = l + aie* r 2 = l + &ie*, (70) 

where \& = (kx + cut) + 9iC and a± and 61 are nonzero even parameters. Introducing T\ 
and r 2 in equation ( 1651) yields the following relation 



W = -a x . (71) 

The dispersion relation 

u + k 3 = 0, (72) 

is obtained from the equation (I64j) . Now, we have ^ = k(x — k 2 £) + #i£ which is similar to 
(|25p but not identical. It is expected since we are interested in travelling wave solutions. 
Since 

U h = F(V) = F(* ) + ^iC^Imo, (73) 

where \l/ = 

t), and u + 6*1^1 = we see that the fermionic solution £1 is 

essentially the derivative with respect to x of u, ie the particular solution of (I69p . In 
the conclusion, we discuss the other solution of this second order linear equation. 

We easily recover a one super soliton choosing a± = i and k — 1. In such a case, 
the Hirota formalism produces the particular travelling wave solution u = (|39|) . 

Now let us exhibit new solutions known as N super soliton solutions [4ll5l [l~2T4T5"] . 

Indeed, in order to find the 2 super soliton solution, we first take 

7i = l + a^* 1 + a 2 e* 2 + ai a 2 A 12 e* 1+ * 2 , (74) 
r 2 = 1 + M* 1 + he"" + M 2 £i2e* 1+ * 2 , (75) 

where now = KjX+a>jt + # 1 £j and and &j are nonzero even parameters (« = 1,2). We 
thus see that we have to take K\ 7^ n 2 , which has the effect of breaking the symmetry 
associated with the translation generator Vt — cV x . Introducing these expressions in 
the mKdV bilinear form ( |64l) and (|65|) . the first equation yields the expected dispersion 
relations: 

Ui + K 3 = 0, i = l,2. (76) 
The second equation gives the following conditions: 

hi = -di, i = 1,2, (77) 

A 12 = B 12 =(^^) 2 (78) 

\«1 + K 2 J 

and a new condition relating the anticommuting variables Ci and C2 given by k 2 7^ 0) 



= «2Cl- 



(79) 



Finally, the r-functions are given by 

n = 1 + a^* 1 + a 2 e* 2 + aia 2 A 12 e* 1+ * 2 , (80) 
t 2 = 1- aie* 1 - c^e* 2 + aia 2J 4 12 e* 1+ * 2 . (81) 

We may enjoy the behavior of the u(x, t) part of the 2-soliton solution 

9(10cosh(|(t-4x)) + 5cosh(t-x) + 8sinh(|(t - Ax)) + 4sinh(t-x)) 



u(x, t) 



72 + 41 cosh(|(3t - Ax)) + 81 cosh(±(7t - Ax)) + 40sinh(|(3t - Ax)) ' 

(82) 



where the parameters k±, k 2 , ai, a 2 are chosen as fti = 2k 2 = 1 and a± = a 2 = z, so that 
= (x — t) + ^iCi and \l/ 2 = |(4x — £) + 0\Q 2 . In Figure 3, we give the solution \u\ for 
t = -10,0,10. 

Again, £ is the derivative of u but we could want to find the other solution of the 
fermionic equation. 



Figure 3. The function \u(x,t)\ for t = -10,0,10. 



For the three super soliton solution, n and r 2 take the explicit form: 

n = l + aie* 1 + a 2 e* 2 + a 3 e* 3 + aia 2 Ai 2 e* 1+ * 2 + aia 3 Ai 3 e* 1+ * 3 

+a 2 a 3 A 23 e M * + a^^is^e* 1 ** 2 ** 3 , (83) 

r 2 = l- a^* 1 - a 2 e* 2 - a 3 e* 3 + a 1 a 2 A 12 e* 1+ * 2 + a^^e* 1 ^ 3 

+a 2 a 3 A 23 e M * - a^a^A^e* 1 ^ 2 ^ 3 , (84) 

where = — + 9iQ (i = 1,2,3) and 

KiQ = Kjd, (86) 

for i, j = 1, 2, 3 (z 7^ j). We may enjoy the behavior of the t) part of the 3-soliton 
solution. In Figure 4, we show again \u\ for t = —15,0, 15. The parameters k±, k 2 , k 3 , 



z. 



ai, a 2 , a 3 are chosen as Ki = yK 2 = ~k 3 = 1 and a\ = a 2 = a 3 

We have thus shown how to construct the N super soliton solution (N = 1,2,3), 
using the bilinear form and fl65l) . by giving the explicit forms of the functions T\ 
and r 2 . The N super soliton solution (iV > 3) is easily generalized using the constraints 
above [17]. We are presently working on a Mathematica program which generates the 
N soliton solution by constructing the r-functions T\ and r 2 for general N. 



Figure 4. The function t)| for t = 



-15,0,15. 



3.2. Rational similarity solutions 

The Hirota formalism can also be used to get rational similarity solutions [10]. We 
assume now a SUSY generalization where the dependent variables T\ and r 2 are 
polynomials in the independent variable z = t~^(x + 9iC) which is similar to ( 1301) . 
To get the SUSY version of the solution (|45|) . we take 

Ti(x,t;9x)=n(z,t) = th, r 2 {x,t;9 1 ) = r 2 (z,t) = t(z 3 + 12), (87) 

since this choice solves the bilinear system (IMj) and (1551) . We thus have the following 
form 

U\x, t;9 1 )= U\z, t) = % log ui - ^ 12 
and, as for the classical case, the invariant solution is thus 

Jb , . „ x . d , /5 3 + 12 



wi(iz) = m x U b (x, t; 6,) = i— log [-^—) ■ (89) 
From a Taylor expansion around #i = we get, = xt~^, 

wx(z) = wtizo) + r^iC-^l^, (90) 

and since -u b = u + 6>i£i = d x U h = t~sw\(z), we get u(x,t) = t~swi(z ) as in ( 1341) and 
£i = C M ^) as expected. 

The generalization to the infinite set of solutions given in [HELL] is direct and we 

where the functions Q n (z) are the Yablonskii-Vorob'ev polynomials define in Appendix 
1. The link with the Hirota formalism is obtained by letting T\ and r 2 be polynomials in 
the independent variable z and the time variable t. In fact, we take the following series 

n(n+l) (n+l)(n+2) 

Tl, n {z,t)=t 6 Q n (2), T2, n (z,t)=t 6 Q B+ i(2) (92) 

which lead to (I9TT) . For example, we see that if n = 1, we recover the expression flHTj) . 



The corresponding invariant solutions w(z) of (I43p are easily obtained by derivation 
,~s ■ d ( Q n+ i(z)\ + 1) 

Wn{z)=t di^\-Qja)> c * = c ^ = *^- ( 93 ) 

where we see that the integration constant c 2 is different for each solution, ie c 2 = c 2;n . 
We also have the following identity 

Ul{z,t) = c 2; „logt + J w n (z)dz. (94) 

We may want to give the corresponding solutions for u n (x, t), in terms of the independent 
variables x and t. For example, from 102(2), we find 

. . 3iUx 2 (x 6 + 24tx 3 + 1440t 2 ) , , 

W2(X ' t] = x* + 72tx* - 8640t3 ' ^ 

and, using w 3 (z), we get 

AiU(x lb + 150te 12 + 9000tV - 604800t 3 a; 6 - 9072000t 4 x 3 - 54432000t 5 ) 



u 3 (x,t) 



x(x 6 + 60tx 3 - 720t 2 )(x 9 + 180tx 6 + 302400i 3 ) 

(96) 



4. Conclusions 

The resolution of the M = 2 SUSY KdV equation for a = —2 has been revised 
with the aim to show some links between the symmetry reduction method and the 
Hirota formalism. In the first case, we have fully used subalgebras and invariants 
including bosonic and fermionic dependent and independent variables to get travelling 
wave solutions and also rational similarity solutions. In the second case, we have been 
able, for the first time, to extend the Hirota formalism to the M = 2 SUSY case and 
to produce solutions of the SUSY KdV equation which are similar to the preceding 
ones but not identical. Indeed, the assumption made has produced a reduction of our 
equation to a M = 1 SUSY one for which the bilinearization is known. From this last 
formalism, we recover the large set of solitonic solutions (N solitons) that has already 
been found in the classical case and in the M = 1 SUSY case. More interestingly, 
similarity rational solutions have been generalized in the SUSY context and we have 
shown that they appear as well when we consider the Hirota formalism. We have shown 
in both cases that £1 = (u x solves the fermionic equation (J69l) for a solution u of 
( I68jl . The Hirota formalism has produced only one independent solution of the second 
order linear differential equation fl69|) . Using classical resolution technics, one can easily 
retrieve the second linearly independent solution of the homogeneous fermionic equation 
for which the Hirota formalism has not yet been adapted to find such solutions. Indeed, 
for i/Q = x — t and u(y ) = sech y , one gets 



6(2/0) = Csech y ( -5 + cosh 2y + 6y tanh y Q ). 



(97) 



For z = xt 3, u(x,t) = t 3Wo(zo) and £1 = t 30 o (z o ) we get 

r,( 2 S 



h ( 6w 2 - \z Q ) (f> = 0. (99) 



which is a solution of the second order linear differential equation, 

1 

3 
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6. Appendix I: On rational similarity solutions of mKdV 

From the approach of Clarkson [TT] adapted to our context, we define the following 
recurrence relation for the construction of the Yablonskii-Vorob'ev polynomials, 

3^Q n+ iQn-i = zQ 2 n - 12 {Q n Q n . zz - Q 2 n . z ) , (100) 

with Qo(z) = 3"^ and Qi(z) = z. Here are some examples, 

Q 2 {z) = z 3 + 12, (101) 
Q 3 ( z ) = 3~s (/ + 60/ - 720) , (102) 

Q 4 ( z ) = - z (z 9 + 180/ + 302400) . (103) 
3 

Using the above polynomials, we can thus construct a series of rational similarity 
solutions of equation (]43p as shown in [TT] , with a slight change of variables: 

±/.a , • d i__ ( Qn{z) 



= (104) 

where the constant of integration C2 is given by 

4,n = Ti^- (105) 

For example, the first solutions w^(z) (w~ = —w+) of the reduce equation (143!) are 
given as: 

= 7, (106) 

jV + 24/ + 1440) 

^2 (*) = 32 „ ■ —7 — — — — , (108) 



2; 



9 + 72/ - 8640 ' 



z 15 + 150/ 2 + 9000/ - 604800/ - 9072000/ - 54432000 
^ 3 W ~ 1 z{z 6 + 60/ - 720) (/ + 180/ + 302400) ' ( ' 



Interestingly, the fermionic equation ( |44p may be solved using some properties of the 
solutions w^. Indeed, let us define 

W ± = -I* + 6(™±) 2 + 6^ 3 (110) 
so that equation (jUJ) with k 2 = becomes 

<« + v£^ = o. (in) 

Depending on the choice of w+ or w~, expression (IllOp takes the form: 

>V+ = jU. + 2 " + 3 + , (112) 

W- = + 12K") 2 + z 2 ^ +3 . (113) 
3 w n + w n+l 

One can see that we have the following property, 

nc +1 = w+, (ii4) 

Z 

for n > and Wq~ = — . We get explicitly, 

3 

z(z 12 + 336z 9 + 2160z 6 + IO8OOOO2 3 + 9849600) . 

WC = — -. (117) 

3 3(2 6 + 60z 3 - 720) 2 1 ; 

From (I114p . we deduce that = 0+ for n > 0. We also get 

<fe + W o -0o = " = 0, (118) 

has a solution expressible in terms of Airy functions. The particular solution <fi =([491. 
solves in fact, 

fc + wr</>r = o. (119) 

7. Appendix II: Bilinear identities 

We present here some identities that are useful to simplify some calculations. In fact, 
taking \&j = KiX + oj{t + 9Q, we get the following identities: 

5P»( e *i. e **) = ( Cl _ C2 + o{ Kl - « 2 ))(«! - K 2 )"e* 1+ * 2 , 
P(X> a! ,A)(e* 1 .l) = 7'(S !B ,fl! t )e* 1 , 



where V is a polynomial. From these identities, we deduce useful results: 

SVf{e^- e* 2 ) = -SVf (e* 2 - e* 1 ), 
SVf +1 (e^- e* 2 ) = SVf +1 (e* 2 - e* 1 ), 
SD^e^-l) = (Ci + ^iXe* 1 , 
SP^e^-e* 1 ) =0. 

For example, if Ti = 1 + <5e* and r 2 = 1 — <5e* with \l/ given as above, we 

(V t + VI){t v r 2 ) = (p t + 1 + 5e*- 1 - SI- e* - 5 2 e*- 

= 25(w + K 3 )e*. 
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